Preprint typeset in JHEP style - HYPER VERSION 



Liverpool Preprint: LTH XXX 



Searching for chiral logs in the static-light decay 
constant. 



UKQCD Collaboration: C. McNeile, and C. Michael 

Theoretical Physics Division, Dept. of Mathematical Sciences, University of Liverpool, 
Liverpool L69 3BX, UK 



Abstract: Using the clover fermion action in unquenched QCD with pion masses as low 
as 420 MeV, we look for evidence for chiral logs in the static-light decay constant. There is 
some evidence for a chiral log term, if the original static theory of Eichten and Hill is used. 
However, the more precise data from the static action of the ALPHA collaboration do 
not show any evidence for non-linear dependence of the static-light decay constant on the 
light quark mass. We make some comments on the connection between chiral perturbation 
theory for decay constants of the pion and static-light meson. 
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1. Introduction 

The data from experiments such as BaBar and Belle is helping to measure the CKM 
matrix better (and hopefully see a breakdown of the standard model formalism). To 
extract information about the quarks, QCD must be solved for various non-perturbative 
matrix elements. In particular, the ratio of the decay constants of the Bg to B mesons 
(^^) is a crucial QCD quantity for the unitarity checks of the CKM matrix. It will become 
more important once Bg mixing is directly measured at run II of the Tevatron. Ali and 
Lubicz Q review the dependence of the QCD matrix elements on the determination 
of I Vtd I and \Vts\. 

There used to a complaisant view (with perhaps a few exceptions [P, ^) that the ratio 
of could easily be computed reliably from lattice QCD, because systematic errors would 
be reduced in ratios of decay constants. 

The error on the ratio of the has recently been increased, however, because the 
uncertainty due to the long extrapolation in the quark mass was underestimated For 
example, the JLQCD collaboration quote Jbs/ fsa = 1-13(3)(1^2)) where the first error 
is statistical and the second error is that from the systematic uncertainties. The dominant 
systematic uncertainty in JLQCD's result is from the chiral extrapolation of /b to light 
quark mass. There has also been work where the ratio of the B meson decay constant to 
the pion decay constant is used to control the log terms ||9|, ^] . 

The problem is extrapolating the value of the /b decay constant from the masses 
in lattice calculations to the physical point. In particular heavy- light chiral perturbation 
theory predicts a log term in the light quark mass dependence of Jb- All previous lattice 
calculations, apart from some preliminary evidence from the HPQCD group 0, |8|, have 
only seen linear dependence of the heavy-light decay constant on the quark mass. 
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There have been a number of attempts to estimate the error from extrapolating down 
to the hght quark masses, using some physically modified form of chiral perturbation 
theory |, |, |, |^]. These have been criticised by Sanz-Cillero et al. who claim that 



the systematic uncertainty due to the chiral extrapolation may have been overestimated by 
JLQCD Q and Kronfeld and Ryan [^. Rather than blindly introducing the chiral log term 
to lattice data when those data do not show any sign of a departure from a linear behaviour, 
it would be better to resolve these issues by explicit looking for non-linear dependence of 
the static-light decay constant on the mass of the light quark and this is the goal of our 
unquenched lattice QCD calculations. 

The UKQCD collaboration have recently finished a calculation that used sea quarks 



with masses around a third of the strange quark mass |12]. This is lighter than the sea 



quark masses used by the JLQCD calculation H . As part of that study UKQCD claimed 



to see some evidence for the chiral log term in the pion decay constant |12|. As has been 
noted by many authors ^] the chiral log structure of f-j^ and Jb is rather similar at one 
loop. Hence, a detection of chiral logs in f.^ is an indication that the parameters of the 
unquenched calculation are close to where chiral logs may occur in the heavy-light decay 
constant. The value of the lightest pion in this calculation is roughly 420 MeV The 



different treatments of the heavy-light chiral perturbation theory of Sanz-Cillero et al. |11 
show that a deviation of linearity is expected at these pion masses. 

The improved staggered formalism has produced gauge configurations with pion masses 
as light as 320 MeV [0]. Heavy-light staggered chiral perturbation theory can produce non- 



intuitive results 1 14]. It is valuable to perform cross checks on results from improved stag- 



gered calculations with another fermion formalism, irrespective of any theoretical concerns 



about the improved staggered formalism |15|. The huge computational costs of unquenched 



calculations with Wilson and domain wall fermions makes this a tough goal |15, 16 1. 
2. Numerical methods 

The basis of our calculation is unquenched gauge configurations generated with the non- 
perturbatively improved clover action and the Wilson gauge action. The lattice parameters 
are: volume 16^ 32, /3 = 5.2, and the clover coefficient was the non-perturbative value of 
2.0171. We only use the same sea and valence quark mass. The full details of the action 
and results on the hadron spectroscopy have been published |12, 17|. 



We use static quarks for the heavy mesons. The static formalism is the ideal framework 
for investigating the log form. It has fewer parameters in the effective Lagrangian (because 



the parameters due to 1/Mq terms are obviously absent). As noted by Arndt and Lin |18 



finite size effects are reduced in the static limit. Also the exploration of chiral logs in 



the light quark sector by CP-PACS |]19[ and qq+q collaboration [^], essentially used a 
relatively coarse lattice spacing. It is very difficult to apply these techniques for heavy 
quark actions because of large aMq errors. 

UKQCD has already published |^l| an extensive analysis of the spectrum of static-light 
mesons and a paper on the mass of the bottom quark |2^. In this paper we look for chiral 
logs in the heavy-light decay constant. 
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As the aim of this work is to look for chiral logs in the /b decay constant that are a 
small effect, it is important to reduce the statistical errors. In our previous calculations we 
were already using all-to-all and fuzzing techniques, so we needed a new method to reduce 
the statistical errors. The number of available gauge configurations is fixed. The ALPHA 
collaboration [^3| have developed a new variant of the static formalism that reduces the 
l/o mass renormalisation that is thought to be the reason for the poor signal to noise ratio 
of static-light calculations. 

The static action is given by 

Sh = a''Y,lp^,{x)DoMx) (2.1) 

X 

DoMx) = -IMx) - W\x - ai)t]Mx - i) (2.2) 
a 

The original static action written down by Eichten and Hill used W{x)t = U{x)t- The 



version of the ALPHA |23| static action that we investigated was: 

3 



V{x)t = - ^[U{x)jU{x + aj)tU{x + at)] + U{x - aj)]uix - aj)tU{x + at - aj)j] (2.3) 
i=i 

where c is a normalisation constant. This is their A variant of the static-light action, here 
labelled by suffix A or called " fuzzed static" . 

The /b decay constant is defined by the matrix element below: 

{0\A^\ B{p)) = ip^fB (2.4) 

The /b matrix element is extracted from the amplitudes in the two point correlator. 

C{t) = Y,{0\A4{x,t)^^j,{x,0)\0) (2.5) 

X 

Z^Zii/g ex]3{—a£t) (2.6) 
where ^'s is the interpolating operator for static-light mesons and ground state dominance 



is shown in equation 2.6. The Zl amplitude is related to the f^^"^ decay constant 



rstatic ry / 2 lystatic lc\ ^7^ 

We discuss the renormalisation factor Z^""^"^^ later. We |^] used all-to-all propagators |24| 
and fuzzed sources to get accurate correlators. We fit a 5 exponential model to a 5 by 5 
smearing matrix. We discuss the fit strategy in more detail in section ^. 

3. Improvement and perturbative matching 

To extract /^*°* we need the renormalisation and 0{a) improvement terms for the static- 
light axial current. 

A'o'"' = i^lkl5^Q (3.1) 



-3- 



The improved static current is: 



(Af"*)o = At' + acrV7fc75^(^fc + ^bV'Q (3.2) 
in the ALPHA formahsm |25]. ^-re covariant derivatives acting on the hght quark 



fields (il)). The improvement term in equation 3.2 was first introduced by Morningstar and 
Shigemitsu p^ . 

To get from the static theory on the lattice to the static theory in the continuum a 
renormalisation factor is required. 

(A|*«*)o = Zl*'^*(l + &i*"*amg)(^f'^*)o (3.3) 

where rriq is the light quark mass. 

The improvement coefficients 6^°* and c^°* have been computed to one loop in per- 



turbation theory |26, 25 1. 



We use the tree level value of 1/2 for ft^''^*"*, because the one loop calculation has not yet 
been completed. 

Ca = -^^-^9 (3.5) 

^A-,stat ^ _Q xig4^2 ^3 g^J 

UKQCD |27] have recently written down the Z factors in the static limit. These were 
obtained from Kurth and Sommer's calculation [pqj. This was an update on the Borrelli 



and Pittori [28| calculation. 



= 1.0 + (i^l^ - 0.137)5^ (3.7) 

As noted by Hernandez and Hill , there is no effect on the value of Za from tadpole 
improvement if the standard exponential fit model is used to extract the amplitude from 
the correlators. In principle the improvement term could be tadpole improved. 

As we are interested in the chiral logs in the leading order heavy-light chiral Lagrangian, 
we don't include a matching factor from the continuum static theory to QCD ||2^, |3^]. The 
MILC collaboration Q have recently discussed the appropriate scale for Z^*"^* using the 
Lepage-Mackenzie scale setting procedure [31|. 

The perturbative expression for Za for the ALPHA static action is not yet available. 
Experience with other "fuzzed" fermion actions suggests that the value of Za will be closer 
to one than for the static heavy quark action of Eichten and Hill ||9|, ^] . The "smearing" 
out of the gauge links can be thought of as averaging over the fields. This reduces the 
perturbative corrections. Alternatively, the "smearing" can be thought of as smoothing 
the potential in which the light quark moves, so reducing the wave-function at the origin, 
related to fs- 
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It was found by the MILC collaboration for a smeared version of the clover 

action, that too much smearing of the fields can drastically change the decay constant. 
For the static action introduced by the ALPHA collaboration, only one level of smearing 
is used. So the heavy quark is restricted to lie within ita of the origin. At our lattice 
spacing, the spatial extent of the light quark in a heavy-light meson is typically 3a (for 
example the node in the excited wave function is at that distance [^4|). Thus we expect this 
smearing of the heavy quark to retain the same qualitative features, but it should affect 
the renormalisation factor for the action. This picture is consistent with the reported 



comparison of fs^ computed in quenched QCD at /3 = 6.2 by Abada et al. ||3^, using both 
the Eichten and Hill static action and the static action introduced by ALPHA. However, 
our results using the Eichten-Hill static quark action should be less contaminated by excited 
states due to the use of all-to-all propagators and variational smearing techniques. 

4. Chiral perturbation theory 

The static limit is the ideal place to study the chiral logs predicted by heavy-light chiral 
perturbation theory. 

'^fB.^fBdVM^, (4.1) 

The one loop correction to the static-light decay constant for 2 flavours of degenerate 
fermions with lightest pseudoscalar mass labelled rriTr is 35] 



?^ = l-^il±^^log(^) (4 2) 



where g is the B*Btt coupling. There have been a number of calculations [36, 37, 33 1 of 
the coupling g using quenched lattice QCD. We use the nominal value of g'^ = 0.35 in this 
analysis. 

In order to compare with our lattice results, we use the expressions derived by Sanz- 
Cillero et al. ||T^ with a fixed cut-off A since this emphasises the region which is reliable in 
chiral perturbation theory. As well as the one loop correction from the lowest order chiral 
Lagrangian, there will be a term from a higher contribution to the chiral Lagrangian, 
namely aim^. Thus there are three free parameters: A, qi and <I>o- In order to establish 
reasonable values for these parameters, we fit this chiral expression at our two heavier 
quark masses, for each of the two values of A = 0.4 and 1.0 GeV, so determining a range 
of predictions for lighter quark masses. 

5. Results 

In table |^ we present our results for the static-light decay constant as a function of sea 
K, value. These values come from fitting our 5x5 matrix or correlators over the t-range 
4-15 with a 5 exponential expression. This sophisticated treatment is optimal to deal with 
the excited state spectrum and to extract cleanly the ground state contribution needed to 
determine Zl. The x^/d-o.f. is acceptable (indeed for k = 0.1350, we also get an acceptable 
for the naive static case for t-range 3-15, and that result is tabulated). 
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Our results show clearly that the fuzzed static source advocated by the ALPHA col- 
laboration does give a much smaller statistical error. 



Name 


No. 


romps 


n 




3/2^ 

^0 Zl 




3/2 A 


DF3 
DF4 
DF6 


160 
119 
139 


1.93(3) 
1.48(3) 
1.06(3) 


0.1350 
0.1355 
0.1358 


0.304lJ^ 
0.282lJ^ 
0.202+^^ 


2.87(11) 
2.89(15) 
2.24(21) 


0.215^^ 
0.179+^ 
0.159+^ 


1.99(3) 
1.79(6) 
1.71(6) 



Table 1: Decay constants (excluding any factors of Za) for the naive static source (Z) and the 
smeared source (Z) for each data set. The uncorrected lattice value is Z^ (in lattice units), while 
the improved decay constant is Zl (in units of ro). 

We compare in fig. |l| the chiral prediction discussed above with our lattice data for 
different quark masses, where the heaviest quark mass corresponds approximately to the 
strange mass. This chiral prediction is arranged to have a slope such that it goes through 
the values at our two larger quark masses for the fuzzed-static case. As can be seen, the 
curvature expected from the chiral logarithm should have set in at our lightest quark mass. 
The statistical errors are not sufficiently small to verify this accurately, although there is 
intriguing sign of a substantial curvature for the naive static case. Since this curvature is 
not reproduced by the relatively more accurate data from the fuzzed static case, we caution 
that it may be a statistical fluctuation. 

If we use the chiral model as a guide to the possible extrapolation to lighter quark 
masses, as discussed above and shown in fig. |l|, we would obtain ^Zl^^b) 1-31 to 1.46 

as A is varied from 0.4 to 1.0 GeV. Thus we have a systematic error of ±5% arising from the 
chiral extrapolation. We also have a statistical error on the slope which is larger at around 
±35%. So we obtain ^ = 1.38 ± 0.13 ± 0.08. This is to be compared with Kronfeld's || 

best guess from JLQCD and HPQCD results of ^ = 1.25 ± 0.10. 

Our primary aim is explore the light quark mass dependence of the heavy-light decay 
constant, rather than to produce new values for the decay constant. As a cross-check on 
our work, we compute the decay constant for static (EH) h quarks. The DF?, data set 
has sea quarks with masses that are close to the mass of the strange quark. As discussed 
in the previous section we use = 0.68. Using tq = 0.55 (0.05) fm, we get /^*"* = 256 
MeV with errors: 4% statistical, 14% from the scale ro and 10% from the uncertainties 
in Za- This value is reasonably consistent with other determinations of /^*"*. ALPHA 
obtain Jb^ = 225 MeV in the continuum limit of quenched QCD using ro = 0.5 fm with 
a non-perturbative renormalisation scheme. Duncan et al. obtain /^*"* = 304 MeV at 



a = 1.78 GeV in quenched QCD using the Wilson action for the light quarks. 

In UKQCD's work on chiral logs in the pion decay constant, there was a concern about 
finite volume effects |12]. At Hgea = 0.1358 it was argued from chiral perturbation theory 
that the finite volume effects were of the order of 8% vo. f-K- A similar order of magnitude 
effect was also estimated by Colangelo and Haefeli [38|. The volume of the lightest data 
;i.5fm)3. 



have studied the effect of the finite 



set DF& is (1.5fm)"^. Recently Arndt and Lin |18, ^ 

volume on the ratio of heavy light decay constants and bag parameters. For a pion mass 
of around 400 MeV in a box of size (1.6fm)'^, they obtain a finite volume effect in the ratio 
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Figure 1: Static-light decay constant as a function of quark mass. 



of 0.006, suggesting that the finite size effects are small. However, the next to leading 
order estimate of finite size effects in /tt was significant |38|. Unfortunately, Colangelo and 
Haefeli [^8| claim that there is not enough information to make a similar estimate for f'B- 

Because of computational limitations we are forced to work at a fixed lattice spacing. 
It is difficult to estimate the uncertainty due to not doing a continuum extrapolation. 
There are now variants of chiral perturbation theory that include the effects of the leading 
lattice spacing errors (see Bar [40| for a review). Aubin and Bernard have developed 
a formalism for heavy-light chiral perturbation theory, at non-zero lattice spacing, with 
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staggered fermions as the light quarks. 



6. Conclusions 

We have looked for the effects of the chiral log in the heavy-light decay constant using the 
lightest unquenched clover sea quarks used to date. Unfortunately, even with pions as light 
as 420 MeV, we do not see any compelling evidence for the chiral log in the static-light 
decay constant. 

We have seen that the quark mass dependence of the static-light decay constant can 
be different to that of the pion decay constant. Both quantities can in principle have a 
different volume dependence and different higher order terms in the chiral Lagrangian. The 
similarity of the one loop expressions for and /b seems to be of limited value 10] in 
guiding extrapolation. 

Although we have been unable to find unambiguous evidence of the chiral logarithm, 
the approach we have tried to take here is potentially a good way to study the light quark 
mass dependence of the heavy-light decay constant. The improvements in the numerical 
techniques for static-light calculations means that quite precise results are achievable. With 
current computer resources, the only hope of working with lighter sea masses with Wilson 
like quarks is to work at relatively coarser lattice spacings. Using the static limit for the 
heavy quarks is more controlled. The use of heavy-light perturbation theory at nonzero 
lattice spacing would also be required. 
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